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ABSTRACT 


Satellite instruments currently provide global maps of surface UV irradiance by 
combining backscattered radiance data with radiative transfer models. The models are 
often limited by uncertainties in physical input parameters of the atmosphere and surface. 
Global mapping of the underwater UV irradiance creates further challenges for the 
models. The uncertainties in physical input parameters become more serious because of 
the presence of absorbing and scattering quantities caused by biological processes within 
the oceans. In this paper we summarize the problems encountered in the assessment of the 
underwater UV irradiance from space-based measurements, and propose approaches to 
resolve the problems. We have developed a radiative transfer scheme for computation of 
the UV irradiance in the atmosphere-ocean system. The scheme makes use of input 
parameters derived from satellite instmments such as TOMS and SeaWiFS. The major 
problem in assessment of the surface UV irradiance is to accurately quantify the effects of 
clouds. Unlike the standard TOMS UV algorithm, we use the cloud fraction products 
available from SeaWiFS and MODIS to calculate instantaneous surface flux at the ocean 
surface. Daily UV doses can be calculated by assuming a model of constant cloudiness 
throughout the day. Both SeaWiFS and MODIS provide some estimates of seawater 
optical properties in the visible. To calculate the underwater UV flux the seawater optical 
properties must be extrapolated down to shorter wavelengths. Currently, the problem of 
accurate extrapolation of visible data down to the UV spectral range is not solved 
completely, and there are few available measurements. The major difficulty is insufficient 
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correlation between photo synthetic and photoprotective pigments of phytoplankton 
absorbing in the visible and UV respectively. We propose to empirically parameterize 
seawater absorption in the UV on a basis of available data sets of bio-optical 
measurements from a variety of ocean waters. Another problem is the lack of reliable data 
on pure seawater absorption in the UV. Laboratory measurements of the UV absorption of 
both pure water and pure seawater are required. 

Keywords: UV irradiance, radiative transfer models, seawater optical properties 

1. INTRODUCTION 

Increased levels of biologically harmful UV-B radiation (280-320nm) resulting from 
the depletion of Earth's ozone layer have been shown to affect aquatic ecosystems. One of 
the important effects of enhanced levels of UVB radiation is a reduction in the 
productivity of phytoplankton caused by inhibition of photosynthesis due to damage to the 
photosynthetic apparatus 1 . Enhanced UVB radiation could also affect the photochemical 
production of carbonyl sulfide in seawater 2 , thereby augmenting the greenhouse effect and 
affecting other long-term global biogeochemical cycles. Photochemical degradation of 
oceanic dissolved organic matter associated with changes in UV radiation flux may affect 
carbon cycling. A detailed overview of the effects of UV radiation on marine ecosystems 
has been published recently 3 . 

The quantitative assessment of UV effects on aquatic organisms on a global scale 
requires an estimate of the in-water radiation field. The total ozone and UV reflectivity 
measurements, from the Total Ozone Mapping Spectrometer (TOMS) satellite 
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instruments, allow calculation of global daily UV irradiance at the ocean surface 4 1 . 
Estimates of UV transmission in ocean waters require knowledge of the inherent and 
apparent optical properties of seawater. For ocean properties, the Coastal Zone Color 
Scanner (CZCS) flown onboard NASA's Nimbus-7 satellite and current ocean-color 
satellite instruments, such as Sea-viewing Wide Field-of-view Sensor (SeaWiFS) and 
Moderate Resolution Imaging Spectroradiometer (MODIS) were designed to provide 
frequent global measurement of water-leaving radiances in the visible region. Seawater 
optical properties and constituents (e.g. chlorophyll concentration) are inferred from the 
water-leaving radiance allowing estimates of inherent optical properties (IOP) in the 
visible region. To calculate the underwater UV irradiance, the visible IOP should be 
extrapolated down to shorter wavelengths. The extrapolation requires some assumptions 
to be justified. 

The main goal of this paper is to assess the problems of the UV penetration into ocean 
waters using global TOMS surface-UV and satellite ocean-color measurements. In 
assimilating these satellite data sets, two major problems arise: the fast radiative transfer 
(RT) modeling of the penetration of UV light into the water and extrapolation of water 
optical properties derived from the satellite visible channels to the UV spectral region. 
The paper discusses both problems. In section 2 we briefly discuss the satellite ocean- 
color sensor and TOMS data as input to the RT models. Section 3 discusses the RT 
models in more details as well as the parameterization of the UV optical properties. 
Section 4 discusses the global products, which can be created from the models. 
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2. SATELLITE DATA 


The Level 3 spatially binned SeaWiFS and MODIS data can be used for estimates of 
chlorophyll concentration and seawater diffuse attenuation coefficient AL(490nm). 
SeaWiFS also provides the daily cloud fraction data. The calibrated radiances (Level- 1) 
over the ocean are atmospherically corrected 8 to derive Level-2 geophysical products, e.g., 
normalized water-leaving radiances, chlorophyll-a 9 , and diffuse attenuation coefficient at 
490 nm, A^490) 10 . These data are spatially binned and averaged on a 9 km global grid 
(Level-3) for each day. 

For underwater irradiance calculations, one needs to know both the direct and diffuse 
components of the surface irradiance, and the boundary conditions at the air-water 
interface. The TOMS standard UV data (described below) provides only the total surface 
irradiance (diffuse plus direct). To calculate the daily-average direct irradiance, 
information on average cloud fraction in each grid-cell is required. Such information can 
be obtained from the 865nm channel of the SeaWiFS sensor. The SeaWiFS 865 nm 
cloud-albedo threshold over the ocean is set at 1.1% albedo". The relation is a binary one: 
if the threshold is crossed, the SeaWiFS pixel is declared cloud contaminated and the 
cloud flag is set for that pixel. It should be noted that the current cloud flag also masks sun 
glint, high aerosols, turbid water, and thin cirrus clouds. A future algorithm will provide 
an improved estimate of the cloud fraction by descriminating the presence of cloud from 
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the presence of heavy aerosol using the 412-490 spectral contrast (described later) 
combined with a technique to examine the spatial variance of the 865 nm reflectance. 

The TOMS daily gridded (Level 3) products (ozone, reflectivity and aerosol index) are 
used as an input to the atmospheric radiative transfer model to generate daily global maps 
of the surface total (direct plus diffuse) spectral irradiance at the satellite overpass time 5 1 . 
To calculate daily UY exposures, diurnal variations of cloud and aerosol amounts are 
neglected. Because of the highly variable nature (temporal and spatial) of cloud cover, the 
TOMS daily UV estimations should be averaged over periods of at least a week to obtain 
a good estimate of the accumulated UV exposure at a specific location 6 . It was shown that 
the corresponded uncertainty in the satellite estimated monthly UV exposure is less than 
5% 12 . 


3. RADIATIVE TRANSFER MODELS 


3.1 Atmospheric model 

The atmospheric RT model provides the boundary conditions at the ocean surface for 
the underwater irradiance calculation. The radiative transfer solutions in the atmosphere 
and in the ocean are coupled through the contribution of photons first reflected from the 
ocean and then scattered back to the water by the atmosphere. However, if the ocean 
albedo is small enough, the atmospheric and oceanic radiative transfer problems can be 
treated separately. The separation of the atmospheric and oceanic RT models gives less 
than 10% resulting error for satellite estimation of underwater UV irradiance 13 . 
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Existing scheme of calculations of surface UV irradiance consists of three steps. The 
first step is calculation of the clear-sky surface irradiance using a lookup table pre- 
computed for pure Rayleigh scattering. Then the clear-sky surface irradiance is corrected 
for non-absorbing aerosols and clouds using a semi-empirical model on the second step. 
The third step is optional; it is performed if absorbing aerosols are detected. Output of the 
scheme is the total (direct plus diffuse) downward surface irradiance 6 . It has been shown 
that the scheme provides reasonable estimates of the total surface irradiance for snow-free 
conditions that compares with ground-based data at 324 nm as well as schemes, which use 
more complicated cloud correction algorithms 7 . However, in the ocean the diffuse and 
direct irradiances are attenuated differently. Therefore, an independent estimation of direct 
and diffuse components is required at the ocean surface. We will briefly describe the 
existing computational scheme and mainly focus on a technique we are proposing to 
estimate the direct and diffuse irradiances separately. 

3.1.1 Clear sky irradiance 

Assuming pure Rayleigh scattering and Lambertian reflection with albedo A at the 
bottom of the atmosphere, the direct and diffuse downward clear-sky irradiance just 
above the ocean surface, Fciear, can be accurately calculated provided the column ozone 
amount is known. In the operational algorithm, Fciear is calculated using Beer’s law for the 
direct component and interpolation from a lookup table of diffuse/direct ratio pre- 
calculated for a Rayleigh atmosphere using climatological TOMS 325DU ozone and 
temperature profiles for different solar zenith angles 6 . Estimates of A can be made from 
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the monthly minimal Lambert equivalent surface reflectivity derived from the Nimbus- 
7/TOMS measurements 14 . For the open ocean regions ^(380nm) typically varies between 
0.05 - 0.08. The satellite measured high-resolution extraterrestrial solar irradiance 
spectrum (the ATLAS-3 SUSIM data) is used in the computations. 

3.1.2 Reduction of UV irradiance by non-absorbing aerosols and clouds 

The common approach for satellite estimations of surface irradiance involves 
calculation of the clear-sky surface irradiance, Fciear , multiplied by Cf. 

Fao, l =Fa«r-C r (1) 

According to the standard semi-empirical model 6 , the factor Ct is a function of the TOMS 
measured scene Lambert Equivalent Reflectivity (LER) at 360 nm, ^ 360 , and surface 
albedo, A , obtained from the minimum LER climatology 14 . This model provides a simple 
algorithm for cloud correction for total irradiance on the ocean surface. To estimate the 
direct and diffuse irradiances separately, we propose using the fractional cloud model 15 , 
with cloud fraction estimated from the SeaWiFS data. The algorithm is as follows. 

First, we estimate the cloud fraction,/ by averaging SeaWiFS cloud fraction data over 
a model grid-cell. For completely cloud-free conditions the TOMS measured LER, i?360 
should be close to the ocean albedo and cloud correction is not required (C r = 1). However, 
due to the possible time differences between TOMS and SeaWiFS overpass (less than an 
hour) and natural geophysical variability in the ocean albedo and cloud amounts, we have 
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to impose a certain threshold on f for clear-sky conditions. Currently, we do not perform 
direct irradiance cloud correction for grid-cells with /<0.05. When f > 0.05, the total 
irradiance is corrected if > A. 

For grid-cells with /> 0.05, an effective cloud reflectivity, R c , is derived from the 
TOMS LER, ocean albedo, A, and cloud fraction,/, using the following expression: 

R = R mZ i WM 
C / (2) 


The Rc is converted to the effective optical depth of the cloud portion of the grid-cell, re- 
using parameterizations based on the radiative transfer calculations . This allows 
calculation of the direct irradiance under the cloud, Fc, direct- 

The grid-averaged direct irradiance is estimated using the following equation: 


P 'direct = fFc, direct + C 1 " /)^0 


(3) 


where Fq is direct irradiance at the surface for a clear sky and Fc, direct is estimated from 
the equation of direct beam attenuation. Finally, the diffuse irradiance is calculated as a 
residue between the total and direct components: 


F, 


Diffuse 


— fr r - F 

” r Clear*' T 1 Direct 


(4) 


where Fciuur is estimated from equation (1), F direct is estimated from equation (3), and Cr is 
estimated from the standard semi-empiric model with replacement of R%o by Rc-- For 



10 


cloud fraction close to 100% the method reduces to the standard TOMS LER method 
with additional direct/diffuse irradiance partition. 

3.1.3 Correction for absorbing aerosols 

An additional correction is needed in the presence of absorbing aerosol plumes, where UV 
irradiance reduction is stronger. The correction is performed using the TOMS aerosol 
index and a semi-empirical conversion factor, which is a function of aerosol height and to 
a less extent of the aerosol type 5 ’ 6 . The absorbing aerosol correction (AAC) algorithm 
accounts for larger attenuation of UV irradiance by absorbing aerosols compare to 
clouds/non-absorbing aerosols of the same reflectivity. For the AAC method the major 
problem arises from uncertainty in aerosol plume height. Current TOMS AAC algorithm 
assumes the nominal height of 3km for plumes of desert dust and biomass burning smoke 
in the tropics 5 ’ 6 . The uncertainty in the actual aerosol height is included in the error budget 
of the TOMS UV product. In the future the aerosol height could be estimated using data 
from assimilation models or other sources. 

While the UV aerosol index from TOMS is the primary tool to correct for the effect of 
absorbing aerosols, the information from a visible channels of SeaWiFS, can provide 
complementary information such as aerosol particle size and optical depth. An algorithm 
to determine aerosol optical depth is already being used operationally with SeaWiFS data. 
An improved version is currently under development that uses the spectral contrast 
between the 412 and 490 nm channels as a discriminator to determine whether an aerosol 
is absorbing or non-absorbing and to choose between different aerosol models 5 . Once an 
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aerosol model has been chosen, radiances from 4 SeaWiFS channels (443, 510, 670, and 
865 nm) are used in a maximum likelihood method to determine aerosol optical depth and 
Angstrom exponent. Results from this algorithm show good agreement with data taken by 
airborne and ground-based sun-photometer measurements during the 2001 Ace-ASIA 
campaign. 

3.2. Radiative transfer in the ocean 

Given the TOMS estimate of the surface UV irradiance, and assuming isotropic 
angular distribution of the diffuse downward radiance at the ocean surface, many 
appropriate radiative transfer schemes can be applied to model light penetration into the 
ocean. There are two basic requirements for those schemes. The RT scheme should be fast 
enough to compute the spectral UV penetration into the ocean on a global scale for 
reasonable time. The RT scheme should have a sufficient accuracy at biologically 
significant optical depths. However, the accuracy of the current optical measurements of 
the fundamental inherent optical properties IOP of seawater is normally about 10%, and 
the errors of extrapolation of these properties into UV spectral region has yet to be 
estimated. The current lack of accuracy in our knowledge of IOP made it reasonable to 
use less sophisticated radiative transfer schemes for the purpose of satellite mapping of 
underwater UV fields 13 . The accurate models 17 ' 18 are very important for validating faster 
algorithms. We are also planning to use the accurate RT code Hydrolight 17 for generation 
of a lookup table of underwater UV irradiance. Interpolation of the lookup table will be 
used in operational algorithms. 
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3.2.1 Fast radiative transfer schemes 

The first model 19 for an assessment of undemater UV radiation and biodoses was 

developed in 1979. In this model it was assumed that irradiance is attenuated 

exponentially with the diffuse attenuation coefficient K d . This simple formulation of the 

radiative transfer in the ocean widely used 20 requires an a-priori knowledge of K d in the 

UV spectral region. In general, cannot be extrapolated from the visible region for use in 

the UV wavelengths. The coefficient K d depends on the angular structure of the light field 

and, thus, on depth (even for a homogeneous ocean), and on seawater inherent optical 

properties (IOPs). Therefore, there is no a-priori reason to expect that K d values in the UV 

region will vary in the same manner with the angular structure of the light field and depth 

as in the visible region. The problem of correlation between spectral values of the diffuse 

*21 

attenuation coefficient has been carefully discussed . 

To calculate UV underwater irradiances, an approximate RT model should have a 
capability to account for the angular structure of the light field. This capability is of 
importance because the direct and diffuse solar fluxes attenuate essentially different. One 
of the approximate RT schemes having this capability is the Quasi-Single Scattering 
Approximation (QSSA) 22 . The QSSA model has a simple analytical formulation, yet 
enabling us to address the dependence of K d on the angular distribution of the light field in 
the ocean. The QSSA is based on the strong absorption with highly anisotropic scattering 
of seawater 22 . It assumes: (a) single scattering in the upward direction; (b) multiple 
scattering in the downward direction in accordance with a delta-function. The spectral 
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irradiance as a function of depth can be written as a sum of the direct solar radiation and 
the integral of surface radiance over spherical angles from the diffuse radiation, both 

n 

attenuated by water ’ . 

The accuracy of the QSSA has been estimated by comparison with the accurate RT 
calculations for the simplified models of the ocean 13,17 ’ 24 . Its accuracy becomes better for 
lower values of the single scattering albedo, normally less than 0.7 in the UV spectral 
region and is even smaller at shorter wavelengths that are more biologically effective. It 
should be also noted that only smaller optical depths play a significant role in biological 
applications of the underwater UV calculations. For optical depths r<5, the QSSA error is 
less than 20% even for high solar zenith angles 24 . All these considerations justified the use 
of the QSSA for calculations of biologically significant parameters from the underwater 
UV irradiance 13 . 

3.2.2. Model of seawater inherent optical properties (IOPs) 

In general RT schemes require the knowledge of all IOPs: the scattering, absorption 
coefficients, and phase scattering function. The QSSA makes use of more limited set of 
IOPs: the absorption coefficient, a, and the backscattering coefficient, bb . The total IOPs 
are the sums of the IOP of pure seawater and scattering and absorbing water constituents: 

a(A.) = a w (X) + a p (X) + a D0M (T); b b (X) = b bw (X) + b bp (A) 

( 5 ) 

where subscripts w, p, and DOM denote the pure seawater, the suspended particulate 
matter (SPM), and dissolved organic mater (DOM), respectively. 
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For a long time, the pure seawater IOPs were usually obtained from a 1981 paper . 
According to recent findings 26 ’ 27 , the pure water absorption coefficient is significantly 
below the previous consensus values 25 in the wavelength range 380 - 500 nm, about 2 
times lower than the old value at 380nm. A recent paper 28 suggests that the most reliable 
combination of absorption data is data 26 for 380 to 700 nm and data 29 for 196 to 320 nm. 
The gap between the data sets, 320 to 380 nm is filled by linear interpolation. It is clear 
from the discussion 28 that additional laboratory measurements and ocean validation are 
needed over the entire UV range. A comparison of the available data sets on the pure 
water absorption in the UV is shown in Fig. 1 . 

The SPM backscattering coefficient and the DOM absorption coefficient can be used 
in the conventional form: 

a Dou W = a o exp[-S(T - . 4, )] ; b hp (T) = 6 0 (/l/T 0 ) _ '” 

( 6 ) 

The DOM spectral slope 5=0.014 nm’ 1 was commonly accepted for the visible spectral 
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region . A more recent study showed that the DOM spectral slope should be made 
slightly greater in the UV spectral region: 5=0.017 ±0.001 nm’ 1 . Recent measurements 
have showed that the DOM spectral slope can increase with photodegradation of colored 
DOM and can vary within a rather wide range from 0.01 to 0.03 nm’ for clear waters . 
Unfortunately, these variations of the DOM spectral slope have not been parameterized. 
Therefore, an average value of the DOM spectral slope for the UV spectral region 
5=0.017 nm' 1 is recommended. The parameter m may vary in a wide range depending on 
the optical type of seawater. Fortunately, the SPM backscattering coefficient, bo, is 
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normally much less than the total absorption coefficient, a, in the UV spectral region. An 

33 

average estimate of the parameter m= 1 is recommended . 

In Case 1 water, where re-suspension of sediments or coastal and terrestrial influences 
are negligible, it has long been recognized that the bulk optical properties are strongly 
correlated with the photosynthetic pigment mass concentrations of the water 33 . The 
quantitative absorption coefficient data combined with photosynthetic pigment mass as 
estimated by chlorophyll-a provide the basis for visible region optical model 
parameterizations 34,35 . The phytoplankton pigment absorption is commonly expressed 
through chlorophyll-a concentration, C, and the chlorophyll-specific absorption 
coefficient: 

a ph (V = Ca' ph (A,C) (7) 

It is well known that the chlorophyll-specific absorption coefficient depends on 
chlorophyll concentration due to, for example, pigment packaging effect. This dependence 
has been parameterized for the visible range 35 : a^ (A, C) = A(A)C R{X \ where the 

functions A(X) and B(A) are tabulated for the visible region. 

The UV particle absorption is more complicated, since there may be strong 
accumulations of pigments with UV absorption 36,37 but only weak correlation with 
chlorophyll concentration. Recent studies of strongly absorbing mycosporine amino acids 
(MAA) indicates that the UV region of the spectrum is not easily modeled based only on 
proxies of bulk photosynthetic pigments, such as chlorophyll-a. Phytoplankton synthesize 
a variety of compounds that absorb radiation in the UV-B and UV-A regions of the 
spectrum and which could affect the response of the cell to UV radiation. In vivo 
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absorption in the UV-A and UV-B shows a wide range of values, with peaks of absorption 
between 320 and 350 nm and a maximum between 330 and 335 nm 38,52 . It is common 
that the in vivo UV absorption is larger than absorption in the blue. The overall variability 
of particulate absorption in the UV is very high compared to absorption in the visible, 
demonstrating that the absorption in the UV is not due to the major photosynthetic 
pigments and that the UV-absorbing compounds, such as MAA, vary independently of 
chlorophyll. This would suggest including an additional term in Eq. 5 that is independent 
of the traditional phytoplankton pigment absorbance. Parameterization of particulate 
matter absorption in the UV and visible has been developed on a basis of data sets of bio- 
optical measurements from the California Cooperative Oceanic Fisheries Investigations 
(CalCOFI). The parameterization follows an approach of Eq. 7. The particulate absorption 
coefficient calculated from the parameterization is shown in Fig. 2 as a function of 
wavelength for different chlorophyll concentration. The correlation between particulate 
matter absorption and chlorophyll concentration is high in the visible but decreases in the 
UV. It has been suggested that UV absorption by MAAs can be estimated using water- 
leaving radiances at 380 and 412 nm 52 . The 380 nm band will be on the GLI sensor to be 
launched on the ADEOS II platform in November 2002. However, the feasibility of 
deriving water-leaving radiances at 380 nm remains to be demonstrated. 

In the recently published paper 13 , we chose a rather simple model of the chlorophyll- 
specific absorption coefficient because of lack of the UV parameterization at that time. 
The model assumes 5(T)= 0 in the UV and adopts the chlorophyll-specific absorption 
coefficient from data 38 . The average for all-stations spectrum was accepted to be the 
chlorophyll-specific absorption coefficient in the model. This parameterization is 
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compared with the CalCOFI parameterization in Fig. 2. The model 13 allows an 
extrapolation of the water absorption and scattering coefficients measured or retrieved 
from satellite measurements in the visible (400 - 600 nm) into the UV spectral region 
(290 — 400 nm). The model contains three input quantities: ao, bo, and C. These 
parameters are to be estimated from available satellite data sets. First, the chlorophyll 
concentration is the standard SeaWiFS and MODIS product. To determine other 
quantities, the Case 1 water model 34 is assumed. According to the model, the DOM 
absorption at 440 nm is 20% of the total absorption of pure seawater and pigments. This 
assumption determines the most important parameter ao . To estimate the backscattering 
coefficient, the standard SeaWiFS product of the diffuse attenuation coefficient and the 
model of the diffuse attenuation coefficient 39 can be used. Estimates of the DOM 
absorption coefficient play the major role in calculations of the UV penetration into 
seawater because backscatter is much less than absorbance in the UV. The largest 
uncertainty is from the way the model of seawater IOPs is constructed. That is, because 
(1) DOM absorption is estimated as 20% of the sum of pure seawater and chlorophyll 
absorption and (2) SPM absorption is calculated using a constant relationship with 
chlorophyll. These two problems should be addressed in further studies. A possible 
approach of determination of DOM absorption can be using empirical algorithms based on 
band ratios 40 . 

Coastal waters are normally referred as Case 2 waters in which IOPs are uncorrelated. 
The above model cannot be directly applied to those waters. However, independent 
retrieval of absorption coefficients of DOM and phytoplankton pigments has been 
suggested 41 ’ 43 . Given the DOM absorption coefficient in the visible region, it can be 
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extrapolated into the UV region. The SPM backscattering coefficient can also be retrieved 
using analytical algorithms 41 ' 43 . 


4. RESULTS 


The biological effect of UV radiation is typically described by action spectra. A large 
number of action spectra, A(k), has been proposed for various biological effects of UV 
radiation in marine environment 44 ' 46 . The biological daily UV doses can be calculated by 
convolution of UV irradiance spectra- Ex(z) with A(X) and integrating over the time of the 
day: 

400 

D(z) = \dt jX (z, e Q {t))A{X)dX (8) 

290 

Comparisons of the simulated dose with measured one were done 13 for the action 
spectrum for unshielded DNA 47 . It was found that calculated daily doses were in a good 
agreement with surface measurements and underwater measurements . 

Using the above-described model, monthly global maps of DNA doses at selected 
depths and 10% penetration depths defined for UVB irradiance and DNA doses were 
created 13 . The main features of the averaged DNA dose map are determined by latitude 
dependence of the surface UV irradiance. The latitude dependence of the DNA dose is 
clearly apparent in all oceans. Some features of the DNA dose map are due to cloudiness 
structure. For example, the cloudiness effect on the DNA dose was observed in the 
Mediterranean Sea - where clear-sky conditions remained for more than a week, resulting 
in DNA dose values characteristic of equatorial regions. It is interesting that the latitudinal 
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distribution of both total ozone and the optical properties of ocean waters are not seen on 
the global DNA dose map. The effects of ozone amount and seawater optical properties 
are almost masked by cloudiness effects, indicating that latitudinal dependence of the UV 
irradiance and cloudiness are the major factors affecting the underwater DNA dose. 
Exceptions to this will be in ocean or coastal areas of large local turbidity. 

An important measure of a biologically weighted dose is the depth at which the dose 
is reduced to 10% of its surface value. This is the approximate depth over which 
biological damage due to UV effects takes place for a particular mechanism. The 10% 
depth depends on the action spectrum used in calculations of UV dose rates. The larger 
the spectral slope of an action spectrum is, the smaller the 10% penetration depth. This is 
because seawater absorbs more strongly in the short-wave region, therefore, shorter 
wavelength radiation penetrates into seawater less than longer wavelength radiation. 

Horizontal distribution of the 10% DNA dose depth is primarily determined by bio- 
optical properties of ocean waters. However, the angular structure of the light incident on 
the sea surface determined by cloudiness structure and solar zenith angle also affects the 
10% DNA dose depth. This is because of the dependence of the diffuse-attenuation 
coefficient on the angular structure of the in-water light field. The 10% UVB irradiance 
depth is normally greater than the 1 0% DNA-dose depth. This is due to the fact that the 
integral over UVB irradiance is mainly determined by the longer wavelength part of UVB 
spectrum as opposed to the DNA dose that is mainly determined by shorter UVB 
wavelengths. The seawater is an effective filter of the shorter UV wavelengths. 

A sensitivity study 13,50 showed that knowledge of the absorption coefficient of pure 
seawater is crucial in estimates of the UV penetration depth. The 10% UVB penetration 
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depth calculated from the extrapolated new absorption coefficients 27 is about 20% greater 
than that calculated from the old coefficients 25 . It is instructive to estimate how variations 
in the DOM absorbance affect the UVB penetration depth. Calculations were conducted 
for two cases. In the first one, no DOM absorption was assumed. The case represents 
upper limit values of the penetration depth. In the second case, it was assumed that the 
DOM absorption at 440 run is 20% of the total absorption of pure seawater and 
pigments 34 . The result demonstrates the significant effect of the DOM absorption on the 
UVB penetration depth. The sensitivity study highlights the importance of accurate 
knowledge of the pure water absorption coefficient and a fraction of DOM absorption. 

5. CONCLUSIONS 

Problems in assessment of the UV penetration into oceanic waters on a global scale 
and some possible solutions were considered. Global mapping of the underwater UV 
irradiance creates challenges for models combining RT computations with assimilation of 
satellite data. The uncertainties in physical input parameters become more serious because 
of the presence of absorbing and scattering quantities affected by biological processes 
within the oceans. We summarized the problems encountered in the assessment of the 
underwater UV irradiance from space-based measurements, and propose approaches to 
resolve the problems. 

We have developed a fast RT scheme for computation of the UV irradiance in the 
atmosphere-ocean system. The scheme makes use of input parameters derived from 
satellite instruments such as TOMS and SeaWiFS or MODIS. The atmospheric part of the 
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model generates spectral direct and diffuse irradiance on the sea surface that are inputs to 
the underwater part of the RT model. The major problem in assessment of the surface UV 
irradiance is to accurately quantify the effects of clouds. Unlike the standard TOMS UV 
algorithm, we use the cloud fraction products available from SeaWiFS and MODIS to 
calculate instantaneous surface irradiance at the ocean surface. Daily UV doses can be 
calculated by assuming a model of constant cloudiness during the entire day. 

The in-water radiative transfer model is based on the QSSA that is simple, 
computationally fast, and yet enables the angular distribution of the light field to be 
addressed. To calculate the underwater UV irradiance the seawater optical properties 
should be extrapolated down to shorter wavelengths. Currently, the problem of accurate 
extrapolation of visible data down to the UV spectral range is not solved completely. The 
major difficulty is insufficient correlation between photosynthetic and photoprotective 
pigments of phytoplankton absorbing in the visible and UV respectively. Empirical 
parameterization of particulate matter absorption in the UV has been done based of 
available CalCOFI data sets. Another problem is the lack of reliable data on pure seawater 
absorption in the UV. Laboratory measurements of the UV absorption of both pure water 
and pure seawater are required. We have developed a simplified model of seawater IOPs 
allowing the extrapolation of the absorption and backscattering coefficients to the UV 
spectral region provided their values in the visible region are known. Values of the 
absorption and backscattering coefficients in the visible region are estimated from the 
SeaWiFS standard products by using the Case 1 water model. 

The sensitivity study has shown that the main parameters controlling levels of the 
most harmful UV-B radiation underwater for clear sky conditions are the solar zenith 
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angle, water bio-optical properties and total ozone. Attenuation of UV-B irradiance and 
DNA dose rate with water depth is primarily controlled by the seawater absorption 
coefficient and its spectral dependence. An influence of the seawater backscatter on the 
attenuation of UV irradiance is considerably less. Changes in the angular distribution of 
the surface radiance due to aerosol load or clouds may result in an irradiance increase (or 
decrease) at a given depth for large solar zenith angles. 

The main spatial features of the monthly maps of underwater DNA dose are 
determined by the SZA and cloudiness. The seawater IOPs and total ozone effects are less 
significant for the spatial distribution of the DNA dose. The spatial distribution of the 10% 
DNA dose depth is mainly determined by the spatial structures of chlorophyll. Cloudiness 
effects and latitude dependence of the 1 0% DNA dose are also observed due to the effect 
of the angular distribution of the light incident on the sea surface in the in-water UV 
irradiance attenuation 
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FIGURE CAPTURES 
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Fig. 1. Comparison of available pure water absorption data sets. Notations: Q&I - data by 
Quickenden and Irvin 29 ; P&F - data by Pope and Fry 26 ; S&F - data by Sogandares and 
Fry 27 ; S&B - data by Smith and Baker 25 . 

Fig. 2. Particulate matter absorption coefficient for different chlorophyll concentrations. 
Solid lines - parameterization based on CalCOFI data, dashed lines - parameterization 
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based on data by Vernet et al. 



ABSORPTION COEFFICIENT, 1/M 


COMPARISON OF PURE WATER ABSORPTION DATA SETS 




ABSORPTION COEFFICIENT, 1 /M 


31 


0.015 
0.010 

0.005 
0,000 

300 320 340 360 380 400 

WAVELENGTH, NM 



0.040 

i— 

q 0.030 
\Z 

Li_ 

0 0.020 
a 

1 0.010 
i— 

CL 

o 0.000 

m 300 320 340 360 380 400 

< WAVELENGTH, NM 



I 0.12 

5 0.10 

C 0.08 

O 0.06 
c_> 

z 0.04 
o 

S: o.o2 

o 0.00 

m 300 320 340 360 380 400 

< WAVELENGTH, NM 


\ 








Problems in assessment of the UV penetration into natural waters from space-based 

measurements 

Alexander P. Vasilkov, Jay Herman, Nickolay A. Krotkov, Mati Kahru, B. Greg Mitchell, and 

Christina Hsu 

Popular Summary 

Increased levels of biologically harmful ultraviolet radiation (UVB, 280-320 nm) resulting 
from the depletion of Earth's ozone layer have been shown to affect aquatic ecosystems. 
Enhanced UVB radiation could also increase the photochemical production of gases 
augmenting the greenhouse effect. Photochemical degradation of oceanic dissolved organic 
matter associated with changes in UV radiation flux may affect carbon cycling. Global 
mapping of the underwater UV radiation creates many new challenges not occurred for global 
mapping of the surface UV radiation. The challenges are mostly related to larger uncertainties 
in physical input parameters caused by biological processes within the oceans. In this paper we 
summarize the problems encountered in the assessment of the underwater UV irradiance from 
space-based measurements, and propose approaches to resolve the problems. 

Satellite instruments, such as TOMS, currently provide global maps of surface UV 
radiation. The major problem in assessment of the surface UV irradiance is to accurately 
quantify the effects of clouds. Unlike the standard TOMS UV algorithm, we propose to use the 
cloud fraction products available from SeaWiFS and MODIS to calculate instantaneous surface 
flux at the ocean surface. To calculate in-water UV radiation we need to know optical 
properties of seawater scattering and absorbing constituents in the UV. Both SeaWiFS and 
MODIS provide some estimates of seawater optical properties in the visible. To calculate the 
underwater UV flux the seawater optical properties must be extrapolated down to shorter 
wavelengths. Currently, the problem of accurate extrapolation of visible data down to the UV 
spectral range is not solved completely, and there are few available measurements. We propose 
to empirically parameterize seawater absorption in the UV on a basis of available data sets of 
bio-optical measurements from a variety of ocean waters. Another problem is the lack of 
reliable data on pure seawater absorption in the UV. Laboratory measurements of the UV 
absorption of both pure water and pure seawater are required 


